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ABSTRACT 

We  present  optical  interferometric  observations  of  the  Be  star  C  Tauri  obtained  using  the  Navy  Prototype 
Optical  Interferometer  (NPOI).  The  multichannel  capability  of  the  NPOI  allows  a  high-quality  internal  calibration 
of  the  squared  visibilities  corresponding  to  the  Ha  emission  from  the  circumstellar  environment.  The  obser¬ 
vations  suggest  a  strong  departure  from  circular  symmetry  and  thus  are  described  by  an  elliptical  Gaussian 
model.  We  use  a  nonlinear  least-squares  fit  to  the  data  to  obtain  the  likeliest  parameters,  and  the  corresponding 
uncertainties  are  determined  using  a  Monte  Carlo  simulation.  We  obtain  3.14  ±  0.21  mas  for  the  angular  size  of 
the  major  axis,  — 62?3  ±  4?4  for  the  position  angle,  and  0.310  ±  0.072  for  the  axial  ratio.  By  comparing  our 
results  with  those  already  in  the  literature,  we  conclude  that  the  model  parameters  describing  the  general 
characteristics  of  the  circumstellar  envelope  of  C  Tau  appear  to  be  stable  on  timescales  of  years.  We  also  compare 
our  results  with  the  known  parameters  describing  the  binary  nature  of  C,  Tau,  and  we  conclude  that  the  envelope 
surrounds  only  the  primary  component  and  is  well  within  its  Roche  lobe. 
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1.  INTRODUCTION 

The  Be  star  C,  Tauri  (=HR  1 9 1 0)  is  one  of  the  brightest  stars 
of  its  type.  This  makes  the  star  a  very  attractive  target  for 
current  optical  interferometers.  In  fact,  the  circumstellar  en¬ 
velope  of  C,  Tau,  with  a  strong  Ha  emission,  has  been  ob¬ 
served  by  several  such  instruments,  including  the  Mark  Ill 
(Quirrenbach  et  al.  1994,  1997)  and  GUT  (Vakili  et  al.  1998). 
This  paper  reports  observations  by  the  Navy  Prototype  Optical 
Interferometer  (NPOI). 

C  Tau  is  a  truly  interesting  target  because  it  shows  very 
significant  variations  in  the  shape  and  the  strength  of  the  Ha 
emission  line.  Although  the  detailed  behavior  of  the  line 
profile  can  be  quite  complex,  the  most  common  spectroscopic 
quantity  used  to  characterize  this  type  of  variability  is  the  ratio 
of  the  violet  to  red  components  {VI R)  of  the  Ha  emission  line. 
This  ratio  can  vary  on  many  timescales,  but  in  the  case  of 
C  Tau,  a  cyclic  VIR  variation  on  timescales  of  a  few  years  is 
observed  (Guo  et  al.  1995).  Most  recently,  this  type  of  vari¬ 
ation  has  been  explained  with  models  having  one-armed 
oscillation  in  a  nearly  Keplerian  disk,  such  as  the  ones  con¬ 
structed  by  Kato  (1983)  and  Okazaki  (1991,  1997),  but  other 
models  have  been  also  proposed  (see  the  models  listed  in 
Guo  et  al.  1995,  and  references  therein).  Although  the  first 
interferometric  studies  related  to  the  variability  of  Be  stars 
suggest  that  one-armed  oscillations  have  been  detected  in  the 
circumstellar  disks  of  C,  Tau  (Vakili  et  al.  1998)  and  7  Cas 
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(Berio  et  al.  1999),  these  first  investigations  do  not  have 
enough  observations  to  model  fully  the  spatial  extent  of  the 
circumstellar  material.  In  fact,  they  rely  on  previously  deter¬ 
mined  model  parameters,  such  as  the  position  angle  (P.A.)  or 
the  axial  ratio  of  the  line-emitting  region,  taken  from  studies 
based  on  observations  from  different  epochs. 

In  this  study  we  model  the  circumstellar  envelope  of  C,  Tau 
as  an  elliptical  Gaussian  and  compare  our  results  with  the 
measurements  reported  by  Quirrenbach  et  al.  (1997)  for  obser¬ 
vations  obtained  from  late  1991  to  late  1992.  We  specifically 
concentrate  on  model  parameters,  such  as  the  axial  ratio,  RA., 
and  angular  size  of  the  major  axis,  that  describe  the  orientation 
and  the  size  of  the  circumstellar  envelope  of  C  Tau.  We  also 
search  for  any  differences  that  might  be  related  to  the  observed 
changes  in  the  shape  and  the  overall  strength  of  the  Ha 
emission  from  different  epochs.  Since  (  Tau  is  also  a  single- 
lined  spectroscopic  binary,  we  also  look  at  how  our  size  of  the 
circumstellar  envelope  compares  with  the  published  binary 
parameters. 

2.  OBSERVATIONS 

The  observations  of  C  Tau  were  obtained  on  1999  March  1 
and  consist  of  eight  individual  scans.  Each  scan  corresponds 
to  90  s  of  integration  during  which  the  fringe  parameters,  the 
squared  visibilities  V^,  are  measured  simultaneously  at  three 
baselines  and  1 1  spectral  channels  covering  the  region  between 
573  and  808  nm.  The  instrumental  configuration  of  NPOI  is 
described  in  detail  by  Armstrong  et  al.  (1998),  and  the  initial 
data  reduction  involved  in  obtaining  the  F  ^-values  from  the 
photon  counts,  for  the  case  of  a  three-beam  combination,  is 
discussed  by  Benson  et  al.  (1997)  and  Hummel  et  al.  (1998).  In 
addition  to  the  observations  of  the  target  star,  we  have  followed 
each  scan  of  (  Tau  with  an  equivalent  scan  on  a  “check”  star 
(for  a  total  of  eight  reference  scans).  The  check  star  was  77  Aur 
(=HR  1641),  which  is  a  B3  V  standard  star  (Garcia  1989)  with 
no  Ha  emission. 

The  circumstellar  envelope  of  a  typical  Be  star  can  emit 
radiation  not  only  in  the  Balmer  series,  but  also  in  other 
spectral  lines,  as  well  as  at  the  continuum  wavelengths  due  to 
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The  entire  calibration  procedure  is  discussed  in  detail  in 
Tycner  et  al.  (2003),  and  we  highlight  only  the  major  steps 
here.  First,  for  a  specific  baseline  and  scan  the  observed 
squared  visibilities  in  the  nine  continuum  channels  are 
expressed  as 

=  (1) 

where  F^  is  a  quadratic  polynomial  representing  the  atmo¬ 
spheric  and  instrumental  degradations,  Vp  is  the  visibility  due 
to  the  stellar  photosphere,  and  Sj  is  the  spatial  frequency  in  the 
radial  direction  from  the  origin  of  the  (M,j;)-plane  for  the  rth 
spectral  channel.  We  assume  that  the  visibility  due  to  the 
stellar  photosphere  is  known  or  can  be  predicted  based  on  the 
estimated  angular  size  of  the  stellar  component.  We  also  as¬ 
sume  that  the  contribution  from  the  circumstellar  envelope  at 
the  continuum  wavelengths  has  a  negligible  effect  on  the 
observed  F^-values. 

The  stellar  photosphere  of  (  Tau  is  nearly  unresolved  with 
our  baselines,  and  therefore  it  is  sufficient  to  model  it  with  a 
simple  uniform-disk  (UD)  component  for  which  the  visibili¬ 
ties  are  given  by 


Fig.  1. — Fourier  («,  o)-plane  coverage  at  Ha  for  observations  of  Tau 
obtained  from  three  simultaneous  baselines  in  eight  scans  on  1999  March  1. 
The  baselines  have  lengths  of  18.9,  22.2,  and  37.5  m,  with  azimuths  of  — 67?5, 
63?6,  and  86?0,  respectively. 


the  hydrogen  free-free  and  bound-free  emission  (Doazan  1 982). 
However,  2  weeks  after  our  interferometric  observations,  the 
Ha  emission  line  of  (  Tau  was  observed  to  have  an  equivalent 
width  (EW)  of  —2.064  nm,  with  a  peak  emission  more  than  3 
times  the  continuum  level  (Banerjee,  Rawat,  &  Janardhan 
2000).  Therefore,  in  our  spectral  region  of  interest,  573- 
808  nm,  we  concentrate  on  the  spectral  channel  that  contains  the 
overwhelmingly  dominant  Ha  emission  (see  also  the  dis¬ 
cussion  in  §  5  of  Tycner  et  al.  2003).  Figure  1  shows  the  Fourier 
(u,  j;)-plane  coverage  corresponding  to  the  Ha  emission  for  all 
scans  used  in  this  analysis.  A  distinctive  pattern  can  be  seen  in 
the  plot  that  is  caused  by  the  changing  projection  of  the  baseline 
vector  on  the  sky  during  the  night  due  to  Earth’s  rotation. 

For  the  observations  presented  in  this  study,  the  alignment 
of  the  spectral  channels  was  such  that  the  Ha  emission  fell 
primarily  within  a  single  spectral  channel.  However,  depend¬ 
ing  on  the  baseline  used,  it  was  either  the  spectral  channel 
with  a  central  wavelength  at  648  or  665  nm  (we  will  use  these 
wavelengths  as  labels  for  the  remainder  of  the  paper).  This 
happened  because,  for  each  baseline,  the  interference  fringes 
produced  at  the  beam  combiner  were  dispersed  onto  a  separate 
lenslet  array,  each  with  slightly  different  spectral  alignment 
(see  also  the  discussion  in  §  3).  The  nine  remaining  spectral 
channels,  both  at  shorter  (four)  and  longer  (five)  wavelengths, 
contain  mostly  continuum  radiation  from  the  star,  and  thus  we 
will  refer  to  them  as  the  continuum  channels. 

3.  CALIBRATION 

The  observed  F^-values  that  we  obtained  after  the  initial 
reductions  are  not  calibrated  with  respect  to  an  external  cali¬ 
brator  and  therefore  are  degraded  by  many  atmospheric  and 
instrumental  effects.  Because  the  observations  are  obtained 
simultaneously  over  all  spectral  channels,  we  can  now  use  the 
continuum  channels  to  calibrate  the  K^-values  in  the  channels 
with  the  Ha  contribution,  which  can  then  be  used  to  model  the 
circumstellar  distribution  of  Ha  emission. 


Vp  =  Fud 


^  Ji  {ttOvds) 

TT01]dS 


(2) 


where  Ji  is  a  first-order  Bessel  function  and  0ud  is  the  angular 
diameter  of  the  UD,  which  for  (j  Tau  is  estimated  to  be  0.39  mas 
based  on  a  photometric  relation  (see  §  3.4  of  Quirrenbach  et  al. 
1997,  and  references  therein).  Therefore,  when  is  known, 
the  coefficients  of  the  quadratic  function  F^  can  be  obtained 
from  the  continuum  channels  using  equation  ( 1 ),  and  then  the 
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Fig.  2. — (a)  Squared  visibilities  of  Tau  obtained  from  11  spectral  chan¬ 
nels  of  the  37.5  m  baseline,  with  significant  atmospheric  and  instrumental 
degradations  still  present.  The  measurements  obtained  from  the  spectral 
channels  at  648  and  665  nm  are  marked  with  open  symbols  and  were  not  used 
in  modeling  the  general  trend  (which  represents  with  a  quadratic 

polynomial,  which  is  shown  as  a  dashed  line,  (b)  Same  data  as  in  (a)  but  after 
the  correction  by  the  T^  factor.  The  predicted  squared  visibilities  for  the 
uniform  stellar  photospheric  disk  are  shown  as  a  solid  line. 
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K^-values  from  the  648  and  665  nm  channels  can  be  calibrated 
according  to  =  F  obs/^^-  In  order  to  analyze  all  three 
baselines  in  the  same  way,  neither  the  648  nor  the  665  nm 
channel  was  used  to  obtain  the  quadratic  forms  of  E^  for  all 
scans,  even  though  only  one  of  these  channels  contains  the 
contribution  from  the  Ha  emission  line  at  a  given  baseline. 

Figure  2a  demonstrates  how  elfectively  a  quadratic  poly¬ 
nomial  {dashed  line)  describes  the  general  trend  seen  in  the 
continuum  channels  {filled  circles)  for  one  of  the  scans  at 
the  37.5  m  baseline.  This  apparent  trend  is  mostly  due  to  the 
atmospheric  and  instrumental  degradations,  but  it  also 
includes  a  slight  contribution  from  the  UD  component.  The 
predicted  F^-values  for  the  UD  component  are  shown  as  a 
solid  line  in  Figure  2b,  and  as  expected,  the  values  of  from 
the  continuum  channels  resemble  those  of  the  UD  model. 

We  can  test  how  well  a  quadratic  function  fitted  to  the 
continuum  channels  describes  the  atmospheric  and  instru¬ 
mental  effects  by  completely  dividing  the  general  trend  (that 
represents  E^  out  of  the  data.  The  normalized  F^-values 
should  cluster  around  unity  within  their  uncertainties,  if  in  fact 
the  quadratic  description  is  satisfactory.  Figure  3  shows  the 
residuals  as  a  function  of  the  central  wavelengths  of  the 


spectral  channels  used  for  all  normalized  observations  at  all 
three  baselines.  When  normalized  by  the  observational 
uncertainties,  the  majority  of  the  residuals  do  fall,  as  expected, 
within  a  band  of  ±  1  ct,  indicated  in  the  plot  by  the  dotted 
lines.  The  only  spectral  channels  that  show  significant  devi¬ 
ations  from  the  general  trend  are  the  ones  with  the  contribu¬ 
tion  from  the  Ha  emission,  the  665  nm  channel  for  the  two 
shortest  baselines,  and  the  648  nm  channel  for  the  longest. 

Figure  4  shows  the  normalized  residuals  for  the  check  star 
rj  Aur  (which  is  also  nearly  unresolved),  and  as  expected,  no 
large  offsets  are  present  in  the  648  and  665  nm  channels.  At 
the  same  time,  small  systematic  offsets  (at  the  level  of  our 
slightly  overestimated  observational  uncertainties)  are  visible 
in  some  spectral  channels.  These  systematic  offsets,  which  are 
spectral  channel  and  baseline  dependent,  also  appear  to  have 
the  same  general  characteristics  for  both  of  our  stars,  and  thus 
one  might  try  to  divide  them  out  from  our  observations  of  the 
science  target.  However,  we  explain  these  systematic  devia¬ 
tions  as  a  limitation  of  how  well  we  can  correct  for  the  ad¬ 
ditive  bias  present  in  each  squared  visibility  determination. 
Based  on  our  experience,  the  additive  bias  can  depend  on  the 
amplitude  of  the  squared  visibility,  and  therefore,  we  decided 
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Fig.  3. — Normalized  residuals  of  Tau  left  after  the  general  trend  present  at  each  scan  and  baseline  has  been  divided  out,  plotted  as  a  function  of  the  central 
wavelengths  of  the  spectral  channels  and  with  the  ±  1  (t  deviation  {dotted  lines).  Each  panel  represents  one  baseline  and  is  labeled  with  the  physical  distance 
between  the  elements  that  define  the  baseline.  The  position  of  the  Ho;  emission  at  656.3  nm  is  indicated  with  a  thick  vertical  line. 
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Fig.  4. — Same  as  Fig.  3,  but  for  the  check  star  r]  Aur 

not  to  divide  these  systematic  offsets  out  of  our  Ha  data  (see 
also  the  discussion  in  §  4.3). 

Because  the  K^-values  due  to  the  UD  component  are  nearly 
the  same  for  the  two  spectral  channels  at  648  and  665  nm,  and 
because  the  channels  have  nearly  the  same  width,  we  can 
combine  the  observations  from  both  to  form  a  single  data  set. 

More  specifically,  we  combine  the  observations  of  (  Tau  from 
the  665  nm  channel  for  the  two  shortest  baselines  with  those 
from  the  648  nm  channel  at  the  longest  baseline,  to  produce  a 
single  data  set  of  F^-values  with  spatial  frequencies 
corresponding  to  the  Ha  line  (i.e.,  we  adopt  a  wavelength  of 
656.3  nm).  Figure  5  shows  the  calibrated  squared  visibilities 
( )  with  the  contribution  from  the  Ha  emission  for  all  three 
baselines  and  eight  scans  (for  a  total  of  24  independent  data 
points),  which  can  now  be  used  in  model  fitting. 

4.  MODELING  THE  Ha -EMITTING  ENVELOPE 
4. 1 .  Best-Fit  Parameters 

The  spectral  channels  at  648  and  665  nm  have  widths  of 
about  16  nm  each  and  thus  are  approximately  7  times  wider 
than  the  Ha  emission  line  (Banerjee  et  al.  2000).  As  a  result,  we 
must  consider  the  contribution  from  the  continuum-emitting 
stellar  photosphere  as  well  as  the  line-emitting  circumstellar 
envelope.  We  treat  this  as  a  linear  superposition  of  the  two 


components,  and  we  model  the  calibrated  -values  shown  in 
Figure  5  as 

VL,.i=[cpVp  +  {l-Cp)K,.]\  (3) 


where  Vp  and  are  the  visibilities  due  to  the  stellar  photo¬ 
sphere  and  the  circumstellar  envelope,  respectively,  and  Cp  is  a 
free  parameter  that  represents  the  fractional  contribution  to  the 
total  flux  of  the  spectral  channel  due  to  the  stellar  photosphere. 

We  model  the  stellar  photospheric  component  in  exactly  the 
same  way  as  we  did  for  the  continuum  channels,  with  a  UD 
that  has  an  angular  diameter  of  0.39  mas  (recall  eq.  [2]).  To 
model  the  circumstellar  envelope  we  use  an  elliptical  Gaussian 
brightness  distribution  centered  on  the  star.  This  is  the  sim¬ 
plest  model  that  resembles  the  power-law  density  distribution 
falloff  usually  assumed  for  the  circumstellar  disks  (such  as  the 
one  introduced  by  Waters  1986)  and  is  also  the  same  model 
used  by  Quirrenbach  et  al.  (1997).  The  normalized  intensity 
on  the  sky  for  an  elliptical  Gaussian  can  be  written  as 


,{0x',dy')  =  ■ 


(Ol/r^  +  e^,) 


1198 


TYCNER  ET  AL. 


Vol.  127 


Fig.  5. — Calibrated  squared  visibilities  of  (  Tau  from  the  spectral  channels 
with  the  Ha  emission  for  all  three  baselines  and  eight  scans.  Shown  are  the 
uniform-disk  model  representing  the  stellar  photospheric  disk  {solid  line),  the 
best-fit  circularly  symmetric  Gaussian  model  {dotted  line),  and  the  best-fit 
elliptical  Gaussian  model  {squares). 

where  r  is  the  axial  ratio,  is  the  angular  size  (FWHM)  of 
the  major  axis  (when  r  <  \),  and  {O^i,  dyi)  are  the  sky  coor¬ 
dinates  in  the  direction  of  the  minor  and  the  major  axes.  The 
transformation  between  these  coordinates  and  the  ones  in  the 
east-west  (9^)  and  north-south  (9y)  directions  on  the  sky  is 


9x  cos  (j)  —  9y  sin  cj), 

(5) 

9x  sin  (j)  +  9y  cos  4>, 

(6) 

where  is  the  RA.  (measured  east  from  north)  of  the  major 
axis. 

The  Fourier  transform  of  equation  (4)  gives  the  visibility 
due  to  the  Gaussian  envelope  and  has  the  form 

Renv(Mx',  v)  =  eXp  +  ^2)/4  ln2],  (7) 

where  Uj^'  and  Vy>  are  the  spatial  frequencies  corresponding  to  9^' 
and  9y',  respectively.  These  frequencies  can  be  easily  expressed 
in  terms  of  the  spatial  frequencies  u  and  v  using  a  coordinate 
transformation  similar  to  that  used  for  the  sky  coordinates  (i.e., 
Ux'  =  u  cos  (j)  —  V  sin  (j)  and  ty/  =  u  sin  (jj  +  v  cos  </>). 

Using  a  nonlinear  least-squares  fit,  based  on  the  Levenberg- 
Marquardt  method  (Press  et  al.  1992),  we  obtained  the  best-fit 
elliptical  Gaussian  model  to  the  data  shown  in  Figure  5  with 
a  best-fit  parameter  set  Rbest  =  (^mj, =  (3.14  mas, 
0.310,  — 62?3, 0.814)  and  a  xl  (per  degree  of  freedom)  of  0.5. 
The  uncertainties  for  these  parameters  are  given  in  the  next 
section  along  with  a  discussion  of  how  they  were  obtained 
using  a  Monte  Carlo  simulation  of  synthetic  data  sets.  We  also 
note  that  the  best-fit  model  parameters  are  influenced  only 
slightly  by  the  adopted  angular  diameter  of  the  UD.  This  is 
because  the  same  model  (Fj,)  is  used  in  the  calibration  of  the 
values  of  F2^(when  the  F^  factor  is  obtained  via  eq.  [1])  and 
in  the  modeling  of  the  Ha -emitting  envelope  (using  eq.  [3]). 


We  have  tested  this  numerically  and  found  that  a  1 00%  change 
in  the  assumed  UD  angular  diameter  produces  changes  in  the 
best-fit  parameters  of  less  than  1%. 

The  elliptical  Gaussian  model  produced  a  much  better  fit  to 
the  data  than  a  two-parameter  circularly  symmetric  Gaussian 
model,  which  resulted  in  an  angular  size  of  1.62  mas  and  Cp  of 
0.69,  but  with  of  5.7.  Figure  5  shows  both  the  symmetric 
and  the  elliptical  Gaussian  models,  along  with  the  observed 
K^-values  for  the  Ha  emission.  It  is  evident  from  the  figure 
that  the  symmetric  model  does  not  match  the  data.  It  is  also 
apparent  that  the  observational  uncertainties  are  slightly 
overestimated,  which  is  also  suggested  by  xl  <  1- 

We  conclude  that  the  elliptical  Gaussian  model  is  a  signif¬ 
icantly  better  match  to  our  data  than  is  the  circular  Gaussian 
model.  The  schematic  representation  of  our  best-fit  elliptical 
Gaussian  model,  as  it  would  appear  on  the  sky,  is  shown  in 
Figure  6.  The  polarization  vector  at  32°  reported  by  McDavid 
(1999)  is  also  indicated  in  the  figure  by  a  dashed  line.  Our 
best-fit  P.A.  of  the  major  axis  is  at  a  right  angle  to  the  polar¬ 
ization  vector,  as  expected  for  the  case  of  electron  scattering  in 
an  axisymmetric  and  optically  thin  disk  (Brown  &  McLean 
1977). 

4.2.  Estimating  the  Uncertainties 

We  used  the  Monte  Carlo  simulation  of  synthetic  data  sets 
(Press  et  al.  1992)  to  obtain  the  estimates  for  the  uncertainties 
of  the  best-fit  model  parameters.  We  started  with  our  original 
data  set  and  obtained  our  best-fit  parameters  for  the  elliptical 
Gaussian  model,  Rbest>  by  minimizing  as  we  have  done  in 
§4.1.  We  then  used  these  parameters  to  create  23,000  syn¬ 
thetic  data  sets,  all  with  the  same  number  of  data  points  as  our 
original  observations,  but  all  with  different  synthetic  measur¬ 
ing  errors.  The  synthetic  errors  were  constructed  by  drawing 
random  numbers  from  normal  distributions,  with  individual 
widths  set  by  the  observational  uncertainties  of  the  original 
data  points.  We  applied  the  same  minimization  scheme  to 
all  synthetic  data  sets  and  obtained  23,000  synthetic  best-fit 
values  for  each  parameter.  This  gives  us  enough  simulated 
parameters  that  the  four-dimensional  probability  distribution 
is  sampled  sufficiently  well,  and  we  are  able  to  measure  the 
different  projected  widths  of  the  distribution  with  high  enough 
precision. 

It  is  not  possible  to  plot  the  distribution  of  all  the  solutions  in 
the  four-dimensional  space,  and  thus  Figure  7  shows  the  mar¬ 
ginalized  distributions  in  six  different  two-dimensional  planes. 
The  majority  of  the  solutions  fall  in  a  localized  region  of  the 
parameter  space,  with  the  highest  concentration  of  solutions 
near  the  region  defined  by  Rbest-  The  density  contours  shown  in 
Figure  7  are  directly  proportional  to  the  probability  density 
functions,  and  therefore  they  provide  us  with  information  not 
only  about  the  widths  of  the  distributions,  but  also  about  the 
degree  of  correlation  present  between  the  individual  pairs  of 
parameters.  We  can  see  from  the  figure  that  most  of  the  pairs  are 
not  correlated,  with  exceptions  in  the  distributions  for  r  versus 
RA.  and  for  Cp  versus  9^j,  where  in  both  cases  there  is  a  de¬ 
tectable  correlation. 

For  the  case  in  which  the  parameters  are  statistically  inde¬ 
pendent  (or  when  the  correlations  between  parameters  are  ig¬ 
nored),  the  two-dimensional  density  contours  shown  in  Figure  7 
can  be  collapsed  further,  by  allowing  only  one  parameter  to 
vary  at  a  time.  This  is  shown  in  Figure  8,  which  displays  the 
individual  probability  distributions  for  all  four  model  parame¬ 
ters.  Because  all  four  histograms  resemble  normal  distributions, 
we  can  fit  a  Gaussian  curve  to  each  distribution  and  obtain  a 
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^  East  6^  (mas)  West  ^ 

Fig.  6. — Schematic  of  the  best-fit  elliptical  Gaussian  model  for  Tau.  The  dashed  line  indicates  the  polarization  vector  at  32°  reported  by  McDavid  (1999).  The 
Roche  radius  of  the  primary  star  (in  the  plane  of  the  disk;  dotted  line)  and  the  projected  circular  orbit  of  the  secondary  component  {dash-dotted  line)  are  both  based 
on  the  orbital  parameters  of  Jarad  ( 1 987).  The  orbit  neglects  the  small  eccentricity  and  assumes  the  same  orientation  of  the  orbital  plane  as  the  circumstellar  disk  (see 
§  5.2  for  more  details). 


corresponding  standard  deviation  (cr)  for  each  distribution. 
Based  on  these  uncertainties,  our  best-fit  elliptical  Gaussian 
model  has  =  3.14  ±  0.21  mas,  r  =  0.310  zE  0.072,  and 
(j)  =  — 62?3  =E  4?4.  We  also  note  that  the  best-fit  fractional 
contribution  Cp  of  0.814zE0.012  compares  well  with  our 
expected  value  of  0.84  zE  0.01,  which  we  estimated  using  the 
EW  of  the  net  Ha  emission  (following  the  same  procedure  as 
Tycner  et  al.  2003)  and  by  approximating  the  stellar  photo- 
spheric  contribution  with  the  NPOl  channel  bandwidth  (i.e.,  by 
adopting  a  stellar  contribution  of  16  ±  0.5  nm). 

4.3.  The  Additive  Bias 

The  data  we  have  presented  in  this  study  might  be  limited 
not  by  the  intrinsic  scatter  of  the  individual  data  points,  but  by 
the  presence  of  an  additive  bias  that  introduces  systematic 
offsets.  The  channel-to-channel  variations  in  our  calibrated 
data  are  at  the  ^1%  level  and  appear  to  be  baseline  dependent 
(recall  Fig.  3).  This  means  that  even  if  our  data  from  the 
spectral  channels  with  the  Ha  emission  are  affected  by  the 
additive  bias,  it  is  unlikely  that  all  three  baselines  are  affected 
in  exactly  the  same  way.  Nevertheless,  to  assess  the  signifi¬ 
cance  of  such  a  potential  systematic  offset  on  our  results,  we 
have  simulated  its  effect  on  our  data. 

We  have  simulated  two  different  cases  where  all  K^-values 
used  in  the  original  model  fitting  (V^^)  were  either  decreased 


or  increased  by  a  constant  offset  of  0.01.  The  best-fit  model 
parameters  were  obtained  from  both  data  sets  using  the  same 
methods  applied  to  the  original  observations,  and  the  results 
are  listed  in  Table  1 .  Although  we  would  intuitively  expect  the 
photospheric  fraction  to  increase  when  the  data  are  shifted  to 
higher  values,  other  parameters  are  also  affected  since  a  linear 
shift  in  F^-values  cannot  be  expressed  as  a  single  adjustment 
of  the  Cp  parameter  (recall  eq.  [3]).  In  fact,  a  constant  positive 
shift  in  K^-values  corresponds  to  larger  changes  at  the  lower 
visibility  values,  thus  making  the  data  look  “less  resolved,” 
and  produces  smaller  values  for  the  angular  extent  of  the 
circumstellar  region,  as  can  be  seen  in  Table  1 .  The  axial  ratio 
and  the  P.A.  of  the  major  axis  are  also  affected  by  the  sys¬ 
tematic  shift,  although  (j)  is  not  very  sensitive  to  the  shift  since 
it  depends  not  only  on  the  F^-values  but  also  on  the  position 
angles  of  the  individual  baselines.  Because  r  depends  on  the 
lower  as  well  as  the  higher  visibility  values  (due  to  the  major 
and  minor  axes),  and  because  the  minor  axis  is  not  signifi¬ 
cantly  resolved,  it  is  affected  slightly  more  by  a  constant  shift 
in  the  data. 

Based  on  the  data  in  Table  1  we  conclude  that  although 
our  results  can  be  affected  by  the  presence  of  an  additive 
bias,  the  magnitude  of  this  effect  is  within  our  quoted  uncer¬ 
tainties.  This  also  justifies  our  use  of  the  overestimated  ob¬ 
servational  uncertainties,  which  yield  <  1  ^nd  give  us  final 
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(a)  (b) 


Fig.  8. — Histograms  of  all  solutions  for  Tau  parameters  plotted  as  a  function  of  only  one  independent  variable.  Each  panel  corresponds  to  23,000  best-fit  values 
obtained  from  the  simulated  data  sets.  Solid  lines  represent  the  Gaussian  fits  to  the  individual  distributions. 

uncertainties  that  also  account  for  the  possible  presence  of  vary  in  any  significant  way  even  if  variability,  such  as  one- 

such  a  systematic  offset.  armed  oscillation,  is  present  in  the  circumstellar  disk.  Our 

best-fit  values  for  r  and  (j)  of  0.310  ±  0.072  and  — 62?3  ±  4?4, 
5.  DISCUSSION  which  we  have  obtained  for  a  single  night,  are  in  excellent 

agreement  with  the  parameters  reported  by  Quirrenbach  et  al. 
5.1.  The  Circumstellar  Envelope  (1997)  for  observations  obtained  over  a  2  month  period  in 

The  circumstellar  envelopes  of  Be  stars  are  typically  as-  1992.  Their  elliptical  Gaussian  model  has  an  axial  ratio  of 

sumed  to  be  thin,  axisymmetric,  circumstellar  disks  viewed  at  0.28  ±  0.02  and  P.A.  of —58°  ±  4°,  as  well  as  an  angular  size 

different  angles,  and  therefore  the  axial  ratio  r  is  mostly  due  of  the  major  axis  of  4.53  ±  0.52  mas. 

to  a  projection  effect  (see,  e.g..  Wood,  Bjorkman,  &  Bjorkman  It  is  not  obvious  that  the  apparent  difference  between  our 

1997).  This  also  means  that  both  r  and  (f)  are  not  expected  to  angular  size  of  3. 14  ±  0.21  mas  for  the  major  axis  and  that  of 


TABLE  1 

Best-Fit  Model  Parameters  at  Different  Levels  of  Simulated  Additive  Bias 


Type  of  Data 

0mj 

(mas) 

r 

0 

(deg) 

Cp 

xl 

. 

3.01  ±  0.24 

0.211 

±  0.120 

-64.5  ±  4.3 

0.817  ±  0.015 

0.5 

f^L . 

3.14  ±  0.21 

0.310 

±  0.072 

-62.3  ±  4.4 

0.814  ±  0.012 

0.5 

eL  -0.01 . 

3.26  ±  0.23 

0.371 

±  0.078 

-60.3  ±  5.2 

0.810  ±  0.012 

0.5 
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Quirrenbach  et  al.  (1997)  is  in  fact  real.  We  have  performed  a 
number  of  tests  to  try  to  assess  the  significance  of  this  ap¬ 
parent  dilference.  We  have  tested  the  null  hypothesis  (that  the 
size  of  the  disk  did  not  change)  by  treating  9^^  as  a  fixed 
parameter  in  our  model-fitting  procedure.  We  used  two  values 
for  0nij  based  on  the  Quirrenbach  et  al.  (1997)  results,  4.53  and 
4.01  mas  (best-fit  value  minus  1  cr),  and  at  the  expense  of 
slightly  worse  fits  (xj  of  1.73  and  1.03,  respectively)  we 
obtained  similar  results  for  r  and  (j)  as  before  (i.e.,  within  their 
quoted  uncertainties).  Although  in  both  cases  the  best-fit 
values  for  Cp  were  slightly  larger,  they  were  in  agreement  with 
our  predicted  value  of  0.84  ±  0.01.  Based  on  this  analysis  we 
conclude  that  we  cannot  confidently  rule  out  the  possibility 
that  the  size  of  the  major  axis  was  the  same  at  the  time  of  our 
observations  and  that  of  Quirrenbach  et  al.  (1997). 

At  the  same  time  there  is  still  a  possibility  (at  the  ~1  a  level 
quoted  by  Quirrenbach  et  al.  1997)  that  this  dilference  is 
caused  by  an  actual  change  in  the  extent  of  the  circumstellar 
disk  or  is  due  to  a  dilferent  spatial  position  of  an  overdensity 
region  such  as  the  one  reported  by  Vakili  et  al.  (1998). 
Although  on  1999  March  16  the  Ha  profile  observed  by 
Banerjee  et  al.  (2000)  does  not  show  a  simple  double  emission 
structure,  it  does  show  clearly  that  the  short-wavelength  side  is 
weaker  than  the  long-wavelength  side  (i.e.,  V jR  <  1).  This  is 
in  the  same  direction  as  the  VIR  of  ~0.7  reported  by  Guo  et  al. 
(1995)  for  the  same  epoch  as  the  observations  of  Quirrenbach 
et  al.  (1997).  Based  on  the  analysis  of  Vakili  et  al.  (1998),  this 
suggests  that  the  overdensity  region  (if  present)  might  be  lo¬ 
cated  in  approximately  the  same  region  of  the  circumstellar 
disk  and  should  affect  the  best-fit  parameters  of  our  simple 
elliptical  models  in  a  very  similar  way.  On  the  other  hand,  it  is 
interesting  to  note  that  the  EW  of  Ha  was  larger  (by  about 
16%)  during  the  observations  by  Quirrenbach  et  al.  (1997) 
than  at  the  time  of  our  observations  (see  the  EW  measure¬ 
ments  by  Guo  et  al.  1995  and  Banerjee  et  al.  2000).  Therefore, 
if  the  apparent  difference  between  our  best-fit  value  for  9^j  and 
that  reported  by  Quirrenbach  et  al.  (1997)  is  real,  then  this 
suggests  that  the  physical  extent  of  the  circumstellar  region 
might  be  directly  correlated  with  the  strength  of  the  Ha 
emission  line.  A  logical  extension  of  our  study  would  involve 
simultaneous  NPOl  and  Ha  spectroscopic  observations. 

5.2.  The  Binary  Component 

C  Tau  is  also  a  single-lined  spectroscopic  binary  with  a 
period  of  about  133  days  and  an  eccentricity  of  ~0.15 
(Harmanec  1984;  Jarad  1987).  The  inclination  of  the  axis  of 
the  orbital  plane  to  the  line  of  sight  is  estimated  to  be 
60°  <  iorb  <  90°,  based  on  the  lack  of  eclipses  and  the  pres¬ 
ence  of  a  hump  in  the  radial  velocity  curve  (see  Floquet  et  al. 
1989  and  references  therein).  A  similar  lower  constraint  can 
be  obtained  on  the  angle  between  the  rotational  axis  of  (  Tau 
and  the  line  of  sight  by  comparing  the  observed  rotational 
velocity  with  the  critical  velocity  predicted  for  a  star  of  the 
same  spectral  type  (Floquet  et  al.  1989).  If  we  assume  that  the 
elliptical  Gaussian  model  we  obtained  in  §  4.1  represents  an 
axisymmetric  disk,  then  the  lower  limit  for  the  inclination 
angle  is  imin  =  arccosr  ~  70°.  If  we  also  assume  that  the 
circumstellar  disk  is  very  thin,  as  suggested  by  Wood  et  al. 
(1997),  we  can  then  use  this  lower  limit  as  an  approximate 
inclination  of  the  circumstellar  disk.  This  inclination  angle  is 
in  the  same  range  as  the  inclination  angles  for  the  star  and  the 
binary  orbit  estimated  by  Floquet  et  al.  (1989)  and  suggests 
that  all  three  planes  have  roughly  the  same  orientation. 


We  do  not  expect  to  detect  any  signature  of  binarity  in  our 
data,  since  the  secondary  component  is  likely  to  be  much 
fainter  than  the  primary.  The  mass  of  the  secondary  is  esti¬ 
mated  to  be  1.0-1. 4  Mq  (Harmanec  1984;  Jarad  1987;  Floquet 
et  al.  1989).  If  it  is  a  1 .4  Mq  main-sequence  star,  it  would  be  an 
F  star  with  My  ^  3  (Gray  1992,  p.  431).  From  the  Hipparcos 
distance  (ESA  1997)  of  128  ±  19  pc  and  the  apparent  V 
magnitude  of  3.0,  the  of  the  primary  is  —2.5  ±  0.3,  or 
5  mag  brighter  than  the  secondary.  For  this  magnitude  dif¬ 
ference,  the  variation  in  as  a  function  of  spatial  frequency 
would  be  at  most  at  the  ^2%  level.  Because  spatial  frequency 
is  proportional  to  2^^  on  a  given  baseline,  such  variations 
would  show  up  as  variations  in  versus  X  in  the  continuum 
channels.  We  see  no  variations  above  the  2%  level  in  any  of 
our  scans,  which  is  consistent  with  the  expected  magnitude 
difference. 

Although  we  cannot  detect  the  secondary  component  of  the 
binary  directly,  we  can  still  obtain  an  interesting  result  if  we 
assume  that  the  inclination  angle  of  the  orbital  plane  of  the 
binary  is  the  same  as  that  of  the  circumstellar  disk  and  use  the 
published  results  from  spectroscopy.  Jarad  (1987)  estimated 
the  basic  parameters  of  the  binary  system  of  C,  Tau  for  various 
choices  of  the  inclination  angle  (assuming  that  Mprimary  ~ 
11.2  Mq)  and  for  i  =  10°  predicts  a  binary  separation  of  254  ± 
20  Rq  and  a  Roche  radius  (RRoche)  of  the  primary  of  144  ± 
12  Rq,  which  is  obtained  when  the  Roche  lobe  is  approxi¬ 
mated  with  a  sphere  of  the  same  volume.  At  the  distance  of 
(  Tau  of  128  ±  19  pc,  the  binary  separation  and  RrocIib  be¬ 
come  9.2  ±  1.5  and  5.3  ±  0.9  mas,  respectively.  Figure  6 
shows  a  schematic  representation  of  the  projected  circular 
orbit  of  the  secondary  component  (we  neglect  the  small  ec¬ 
centricity)  and  the  extent  of  the  Roche  radius  in  the  plane  of 
the  disk,  along  with  our  best-fit  elliptical  Gaussian  model. 

The  inferred  binary  separation  is  large  enough  to  place  the 
secondary  component  outside  the  circumstellar  disk,  in 
agreement  with  many  models  (see,  e.g.,  Castle  1977).  Fur¬ 
thermore,  the  disk  lies  well  within  the  Roche  radius,  since  the 
outermost  disk  radius,  which  we  can  take  to  be  the  1%  level  of 
our  Gaussian  model,  is  ~4  mas,  which  suggests  that  the  disk 
may  be  truncated  by  tidal  interaction  with  the  secondary.  If 
this  is  the  case,  we  would  not  expect  to  find  any  strong  radio 
emission  from  a  region  much  more  extended  than  the  Ha- 
emitting  region  (in  fact,  Taylor  et  al.  1990  did  not  detect  any 
radio  emission  for  the  (  Tau  system).  The  (  Tau  system  might 
be  one  of  the  few  Be  stars  for  which  the  presence  of  a  nearby 
companion  limits  the  extent  of  the  circumstellar  disk. 

6.  SUMMARY 

We  have  fitted  an  elliptical  Gaussian  model  to  interfero¬ 
metric  observations  of  C  Tau  obtained  on  1999  March  1  with 
the  NPOI  and  have  shown  that  it  effectively  describes  our 
data.  We  have  also  shown  that  a  circularly  symmetric 
Gaussian  model  is  not  capable  of  modeling  the  data.  Our  best- 
fit  values  for  the  axial  ratio  and  RA.  of  the  elliptical  Gaussian 
confirm  the  values  reported  by  Quirrenbach  et  al.  (1997)  for 
observations  obtained  about  7  years  earlier.  We  conclude  that 
the  orientation  of  the  disk  does  not  appear  to  change  on 
timescales  much  longer  than  the  orbital  period  of  the  binary. 
We  also  conclude  that  the  circumstellar  envelope  is  smaller 
than  the  size  derived  by  Jarad  (1987)  for  the  Roche  lobe  of  the 
primary  component  of  the  binary. 

Although  our  best-fit  value  for  the  major-axis  angular  size 
appears  to  be  smaller  than  the  value  reported  by  Quirrenbach 
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et  al.  (1997),  obtained  when  the  Ha  emission  was  stronger,  the 
difference  is  only  slightly  larger  than  the  sum  of  the  formal 
errors,  and  therefore  it  is  difficult  to  assess  the  significance  of 
this  change.  At  the  same  time,  if  this  difference  in  9^^  is  real,  it 
suggests  that  the  size  of  the  circumstellar  envelope  is  correlated 
with  the  strength  of  the  Ha  emission.  Further  interferometric 
observations  of  this  system,  preferably  using  the  same  instru¬ 
mental  setup,  should  be  capable  of  investigating  the  relation¬ 
ship  between  the  Ha  emission  strength  and  the  characteristics 
of  the  circumstellar  envelope. 

The  variability  of  ^  Tau  is  known  from  spectroscopy  (see, 
e.g.,  Hanuschik  et  al.  1996),  as  well  as  from  optical  interfer¬ 
ometry  (Vakili  et  al.  1998),  and  thus  it  is  important  to  acquire 
sufficient  data  to  model  the  source  characteristics  on  short 
enough  timescales.  We  have  demonstrated  in  this  paper  that 
satisfactory  model  fits  can  be  obtained  from  only  12  minutes 
of  data  using  a  three-beam  configuration,  and  therefore,  the 
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NPOI  is  capable  of  detecting  intensity  variations  in  the  cir¬ 
cumstellar  environments  of  Be  stars  on  short  timescales.  This 
is  an  important  and  unique  observational  tool,  especially  if 
combined  with  simultaneous  high  spectral  resolution  Ha 
observations,  that  can  probe  new  astrophysical  phenomena. 
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